Long-range antiferromagnetic interactions in ZnFe 2 04 and CdFe 2 04 
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For the first time, the Fe-Fe interactions in the geometrically frustrated antiferromagnetic systems 
of zinc and cadmium ferrites are determined quantitatively by the first-principles methods of density 
functional theory. Both the generalized gradient approximation (GGA) as well as GGA plus the one- 
site Coulomb interaction (GGA+U) are considered for the exchange-correlation energy functional. 
The interactions up to third neighbours are found to be all antiferromagnetic for both materials, 
regardless of which approximation scheme ( GGA or GGA+U) is used. Surprisingly, the third- 
neighbour interactions are estimated to be much stronger than the second-neighbour interactions 
and on the same order in magnitude as the first-neighbour interactions. 

PACS numbers: 61.50.-f,71.20.-b,75.50.Ee 



The geometrically frustrated antiferromagnets, dis- 
tinct from other magnetic systems, exhibit a number of 
exceptional features whose formation, to a great extent, 
depends on the range as well as the nature of the interac- 
tions between the magnetic ions in the systems |l|, [2|, |3j. 
Both zinc and cadmium ferrites (ZnFe2 04 and CdFe2C>4 
shorted as ZFO and CFO hereafter) are well known as 
one type of the geometrically frustrated systems 0, Q. 
Magnetic neutron scattering measurements have been 
carried out on both ZFO and CFO to determine the in- 
teractions between the Fe ions. Surprisingly, it was con- 
cluded that the frustration in CFO [6] is driven mainly 
by the strong nearest-neighbour antiferromagnetic inter- 
action in contrast to the antiferromagnetically coupled 
third- neighbour interaction in ZFO[7|. In this Letter, we 
study the magnetic interactions between the Fe ions in 
these two compounds by considering structures with dif- 
ferent collinear magnetic distributions as well as cation 
distributions using ab initio methods. 

The structure of spinel ferrites are constructed by fill- 
ing one-eighth of the tetrahedral sites and half of the 
octahedral sites (denoted as A and B site respectively 




FIG. 1: The four B-site cation distributions considered in 
this study. The numbers denote the nth layer when one 
views along the (001) direction of the cubic cell. The B-site 
cations are divided into two groups (Bl and B2 in Tab. I), 
one denoted in green (light grey) and the other in blue (dark 
grey) . The two types of B-site cations correspond to different 
cations, e.g. Fe and Zn, in the inverse spinel or oppositely 
polarized Fe ions in the normal spinel. 



hereafter) in the FCC sublattice of oxygen. Spinel fer- 
rites are usually categorized into two types, i.e. normal 
and inverse spinel, according to the distribution of the 
divalent and trivalent cations in the A and B sites. A 
normal spinel corresponds to the structure with all the A 
sites being occupied by the divalent cations, e.g. Zn (Cd) 
of ZFO (CFO), and the B sites by the trivalent cations, 
i.e. the Fe ions. In an inverse spinel structure, the di- 
valent cations occupy half of the B sites while the other 
half as well as all the A sites are occupied by the Fe ions. 
Both ZFO and CFO are considered forming normal spinel 
though coexistence of both phases has been reported Q]. 
To investigate the magnetic interactions between the Fe 
ions, we consider three different B-site cation distribu- 
tions (named SCI, SC2, and SC3 as shown in Fig. 1) 
on the basis of a cubic cell consisting of eight formulas 
(56 atoms) of the materials. Inclusion of the possible 
cation distributions on the A and B sites as well as the 
different magnetic orders of the Fe ions leads to the struc- 
tures summarized in Tab. I, i.e. four configurations for 
normal spinel and 6 configurations for inverse spinel. A 
fourth B-site cation distribution, i.e. SC4, is also in- 
cluded to estimate the magnitude of the more distant in- 
teractions between Fe cations in the normal spinel as dis- 
cussed later. There have been first-principles calculations 
for other spinel compounds, e.g. studies of MnFe204, 
Fe 3 4 , CoFe 2 4 , and NiFe 2 4 by Szotek et alQ and 
ACr 2 X 4 (A=Zn, Cd, Hg; X=0, S, Se) by Yaresko^. 
Note that the o nly configuration considered in the pre- 
vious studies [J, [ll[ corresponds to that of nbl in Tab. I, 
i.e. the normal spinel phase with an equal number of 
oppositely polarized B-site Fe ions distributed as those 
in SCI of Fig. 1. We shall show that the nb3 is a more 
stable structure compared to that of nbl and similarly 
ib3 to ibl for both ZFO and CFO. 

To study the extent of the interactions between the 
cations in these systems, the distances and numbers of 
neighbouring cations for the A- and B-site cations with 
neighbouring distances up to around 6.7 A are listed in 
Tab. II. If we assume that the interactions between the Fe 
ions in the normal spinel can be approximately modelled 



TABLE I: The ten configurations considered in the present 
study and their corresponding cation distributions on the A 
and B sites as well as the magnetic orders of the Fe ions. 





site 


A Bl B2 


configurations 


normal 


cation 


Zn(Cd) Fe Fc 


spinel 


magnetic 
configurations 


+ + 
+ 


nal=na2=na3 
nbl,nb2,nb3 


inverse 


cation 


Fe Zn(Cd) Fc 




spinel 


magnetic 
configurations 


+ + 
+ - 


ial,ia2,ia3 
ibl,ib2,ib3 



TABLE II: The distances (d) and numbers (N) of neighbour- 
ing cations in the spinel systems are listed in the order of 
increasing distances. For example, the number of the nearest- 
neighbour B-site cations for an A-site cation is twelves as in- 
dicated in the second line. The distances are written in terms 
of the squared values (in unit of ao, i.e. the cubic lattice con- 
stant) and the length (in unit of A) assuming ao = 8.5 A. They 
are further grouped according to the order of neighbouring. 
The nth-neighbour interactions between the B-site Fe ions in 
normal spinel, i.e. J n of the Ising model, are also denoted. 
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by the Ising model, i.e. H = -(1/2) J2nJ2i Jn,SiS i+n , 
the values of Eq, Ji, J2, and J3 of these materials can 
be obtained from the calculated total energies of the four 
configurations of na, nbl, nb2, and nb3. The approxi- 
mations made by using this simple model to describe the 
interactions in these compounds and to obtain the corre- 
sponding interaction energies are supported by the fact 
that the moments on the Fe ions in ZFO and CFO are 
sizable and considerably localized. Similar approaches 
have been used previously to study the magnetism in, e.g. 
spinel MnC^PJj and YCu02.5[13|]. Studies using first- 
principles calculations for these systems with collincar 
mangetic distributions should provide information on the 
size and the extent of the interactions between the mag- 
netic cations in these systems. These information are, on 
the other hand, important references to further under- 
stand the possible phases of the materials through the 
studies of the corresponding classical three-dimensional 
Heisenberg model. 

All electronic calculations in this study are based on 
the spin-polarized density functional theory [lj, [l5| . The 
generalized gradient approximation (GGA) proposed by 
Perdew, Burke, and Ernzerhof 16] for the non-local cor- 



rection to a purely local treatment of the exchange- 
correlation energy functional are used. The on-site 
Coulomb interaction U ,17] for Fe ions is also included 
in the second stage of the investigation to study its ef- 
fect on the physical properties of these systems. The 
values of U and J are taken from the previous studies as 
4.5eV[ll] and 0.89eV[I3] respectively. The interactions 
between the ions and valence electrons are described by 
the projector augmented-wave (PAW) method [19] in the 
implementation of Kresse and Joubert [20| and the num- 
bers of the treated valence electrons are 12, 8 and 6 for Zn 
(Cd), Fe and O atoms respectively. The single-particle 
Kohn-Sham equations [21| are solved using the plane- 
wave-based Vienna ab-initio simulation program (VASP) 
developed at the Institut fur Material Physik of the Uni- 
versitat Wien [22| . The energy cutoffs used for the plane- 
wave basis is 500 eV and a set of (6 6 6) k-points sampling 
according to Monkhorst-Pack [23| is used for the integra- 
tion over the first Brillouin zone, unless specified other- 
wise. Relaxation processes in optimizing static structures 
are accomplished by moving oxygen atoms to the posi- 
tions at which the atomic forces are smaller than 0.02 
eV/A. The cell volumes are also relaxed under the con- 
straint that the systems remain cubic. 

The calculated energies (relative to the corresponding 
lowest-energy one), lattice constants and magnetic mo- 
ments of Fe ions of the studied configurations (Tab. I) 
for both ZFO and CFO using either GGA or GGA+U 
approach are listed in Tab. III. 

The experimental lattice constants for ZFO and CFO 
are around 8.52AQ and 8.72A@ respectively. The cal- 
culated lattice constants are all within 1% deviation from 
the experimental values except for the ial, ia2, and ia3 
of ZFO [24]. One also notices that the energies of these 
three ferromagnetic inverse-spinel configurations are sig- 
nificantly larger than the rest of the considered configu- 
rations. Therefore the effect of U is not applied to the 
ial, ia2, and ia3 configurations. From now on we shall 
focus the discussions on the rest seven configurations. 

The local magnetic moments for Fe ions in these config- 
urations were found to be within the range of 3.7~3.9/io 
and 4.1^4.2 ^ for GGA and GGA+U calculations re- 
spectively. The GGA+U results are closer to the exper- 
imental estimates of the magnetic moments, i.e. 4.22|25| 
and 4.440/l/o for ZFO and CFO respectively. 

The orders in energy of these configurations are differ- 
ent for ZFO and CFO which demonstrates the underly- 
ing differences in the properties of these two compounds. 
However the lowest energy configurations are identical for 
these two compounds, i.e. they are ib3 and nb3 in the 
GGA and GGA+U calculations respectively. Alhough 
the lowest energy configuration is an inverse spinel in the 
GGA calculations, it is a normal spinel in the GGA+U 
calculations. The on-site Coulomb interaction U is usu- 
ally considered as essential for properly describing the 
electronic and magnetic properties of the Fe ions in spinel 
ferrites. Note that the lowest energy configuration of the 
GGA+U results, i.e. nb3, is different from the configu- 



TABLE III: The calculated energies (in unit of eV per 8 formulas of the materials) relative to the corresponding lowest-energy 
configuration, the lattice constants (in A) and the local magnetic moments of Fe ions (in unit of Bohr magneton) for the ten 
configurations of Tab. I using either GGA or GGA+U (data in parenthese) approach. 



ZFO 


na 


nbl 


nb2 


nb3 


ial 


ia2 


ia3 


ibl 


ib2 


ib3 


AE 

Oo 


4.28(0.54) 

8.52(8.53) 

3.9(4.2) 


2.21(0.18) 

8.50(8.52) 

3.8(4.2) 


2.71(0.14) 

8.51(8.52) 

3.9(4.2) 


1.98(0.00) 

8.50(8.52) 

3.8(4.2) 


5.18 

8.35 

3.7/1.5 


4.89 

8.36 

3.7/1.4 


4.12 

8.35 

3.7/1.5 


0.38(2.18) 

8.45(8.49) 

3.7(4.1) 


0.90(2.68) 

8.45(8.50) 

3.7(4.1) 


0.00(1.68) 

8.45(8.49) 

3.7(4.1) 


CFO 


na 


nbl 


nb2 


nb3 


ial 


ia2 


ia3 


ibl 


ib2 


ib3 


AE 

Oo 


3.50(0.76) 

8.82(8.82) 

4.0(4.2) 


0.99(0.15) 

8.78(8.81) 

3.8(4.2) 


1.74(0.27) 

8.80(8.81) 

3.9(4.2) 


0.79(0.00) 

8.79(8.81) 

3.8(4.2) 


8.32 

8.77 
3.8 


7.71 
8.15 
3.9 


6.89 
8.77 
3.8 


1.91(4.94) 

8.73(8.76) 

3.7(4.1) 


1.89(4.41) 

8.75(8.78) 

3.7(4.1) 


0.00(2.81) 

8.72(8.75) 

3.7(4.1) 



TABLE IV: The nth-neighbour interaction energy (in unit of 
meV) between the B-site Fe ions obtained from the calculated 
energies of the normal spinel configurations of na, nbl, nb2 
and nb3. The J\ and J 2 obtained from considering merely the 
energies of nbl, nb2 and nb3 are also listed. Unless specified 
accordingly, the calculations use an energy cutoff of 500eV 
and the Monkhorst-Pack mesh of (6 6 6) for the k-points 
sampling. 



ZFO 


4Ji 


4J 2 


4J 3 


4Ji Q 


4J 2 b 


PBE 

PBE+U C 


-220.4 
-97.7 


-19.0 
-7.0 


-23.6 
-16.5 


-126.2 
-31.6 


+28.1 
+26.0 


PBE+U" 

PBE+U 

PBE+U e 


-43.7 
-39.6 
-39.7 


-3.2 

-2.7 
-2.6 


-12.2 
-12.3 
-12.4 


+5.0 
+9.7 
+9.8 


+21.1 
+22.0 
+22.2 


PBE+U 1 
PBE+U ff 


-38.1 
-40.2 


-3.5 

-2.4 


-12.7 
-12.5 






CFO 


4Ji 


4J 2 


4J 3 


4Ji a 


4J 2 b 


PBE 
PBE+U 


-275.9 
-71.1 


-18.7 
-2.2 


-21.8 
-10.6 


-188.8 
-28.9 


+24.8 
+18.9 



a ,h = (nbl - nb2)/2 b J 2 = (nbl - nb3)/4 c U=2.5eV 
d E cuto// =400eV kpoint set=(4 4 4) 

e E cuto //=600eV kpoint set=(8 8 8) 
* obtained from na, nbl, nb2, nb4 
9 obtained from na, nbl, nb3, nb4 



ration considered in all the previous studies, i.e. nbl. 

The interaction energies between the B-site Fe ions up 
to third neighbours can be determined by applying the 
Ising model to the calculated energies of the four normal 
spinel configurations, i.e. nbl, nb2, nb3 and na, and the 
results are presented in Tab. IV. They all turn out to be 
antiferromagnetic, no matter whether the effect of U is 
included or not. The effect of U reduces the magnitudes 
of J\ and J2 substantially, but only halves the magni- 
tudes of J3. The universal reduction in magnitude of J n 
after including U is implied in the density of states of 
these systems. The most prominent effect of U in ZFO 
is localizing the five d electrons of Fe ions which are pre- 
viously well hybridized with the other valence electrons, 
and therefore weakening the interactions between them 
while at the same time widening the band gap of the sys- 
tem. Similar results take place in CFO except that the 



strongly localized ten d electrons of the Cd ions locate 
at the lower energy range than the localized five d elec- 
trons of Fe ions even after the effect of U is included. 
That how the exchange interactions depend on the val- 
ues of U have been studied previously, e.g. in FeSb0426] 
and in NiGa2 84(23. The variations of J due to U were 
shown being monotonically. Our results of using U=0 
and U=4.5 suggest that for both ZFO and CFO with 
the on-site Coulomb interaction within the ragne of to 
4.5 all lead to having antiferromagnetic interactions up 
to third neighbours and stronger third-neighbour interac- 
tions than the second-neighbour ones. These conclusions 
are consistent with the tested results of calculations using 
U=2.5eV as shown in Tab.IV. 

As the magnitude of J2 in the GGA+U results is con- 
siderably smaller than those of J\ and J3 , the numerical 
errors due to using finite basis sets and discrete k points 
are examined. The lower (400eV) and higher (600eV) 
energy cutoffs for the plane-wave basis are considered as 
well as the sparser (4 4 4) and denser (8 8 8) sets of 
Monkhorst-Pack sampling. The results in Tab. IV in- 
dicate that the significant figure for the values of these 
interactions, i.e. 8 J n (8 is the number fomulas in the cu- 
bic cell used in the calculations), can be considered as in 
meV. The magnitude of J3 is at least three (five) times 
larger than that of Ji and about one third (sixth) of J\ 
for ZFO (CFO) in the GGA+U calculations. This, i.e. 
the J3 interaction, corresponds to an interaction between 
B-site Fe ions of at least 6 A apart which is a considerably 
long-range interaction. If the interactions beyond the 2nd 
neighbours are negligible, the values for J\ and J2 can be 
readily obtained from merely considering the energies of 
nbl, nb2 and nb3. These results are also listed in Tab. 
IV. Under this assumption, both Ji and J2 of ZFO in 
the GGA+U calculations are switched to ferromagnetic 
interactions. This is not a plausible result when con- 
sidering the experimentally identified antiferromagnetic 
phase for ZFO. To estimate the magnitude of the next 
distant interaction, i.e. J4, which is not included in the 
previous discussion, the total energy of the SC4 configu- 
ration in the normal spinel, i.e. nb4, was also calculated. 
The magnitude of J4 can not be obtained from consider- 
ing the five configurations of na, nbl, nb2, nb3 and nb4 
as their energy fomulas are linearly dependent. How- 



ever, the magnitudes of the first three neighbour interac- 
tions can be evaluated from considering cither na, nbl, 
nb2 and nb4 or na, nbl, nb3 and nb4 whose results are 
also listed in Tab. IV. It is obvious that the effect of the 
longer-distance interactions does not change the findings 
we shall conclude next, i.e. existing antiferromagnetic 
interactions up to the third neighbours, much weaker Ji 
compared to those of J\ and J3 and that J\ and J3 are 
at the same order of magnitude. Similar characters take 
place in CFO except that the nearest-neighbour interac- 
tion, when compared to J3, is much more dominant than 
that in ZFO which is consistent with the experimental 
suggestions [a 0| . 

There have been numerous theoretical studies for the 
geometrically frustrated antiferromagnetic systems using 
the classical three-dimensional Heisenberg modelp], 0, y]- 
Studies including up to the 2nd neighbour interactions 
suggested that with the present estimated values of Ji 
and J 2 , i.e. J\ < and J2/J1 < 0.5, the system re- 
mains paramagnetic down to 0K[2J. A mean-field ap- 
proach to magnetic ordering in the highly frustrated py- 
rochlores also concluded that no long-range order can be 
established for a system with antiferromagnetic J\ and 



J3 interactions but negligible J2 and J4 interactions Q. 
Experimentally there has been suggestion that no long- 
range order can be established down to 0.1K for CFO 
even with an applied field of up to 9T[6| and similarly 
ZFO was found remaining disordered even at the lowest 
observable temperature of 1.5K 7]. 

In summary, we have generated different configura- 
tions, in both cation distribution and magnetic order, to 
study quantitatively the Fe-Fe interactions in ZFO and 
CFO using GGA as well as GGA+U approach. The in- 
teractions bewteen the B-site Fe ions are found to extend 
up to 3rd neighbours and all in antiferromagnetic nature. 
The 2nd-ncighbour interaction is estimated to be much 
smaller in magnitude than those of the first-neighbour 
and the third-neighbour interactions. 
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